arXiv: 1503.0247lvl [quant-ph] 9 Mar 2015 


Geometry of Gaussian quantum states 


Valentin Link and Walter T. Strunz 
Institut fur Theoretische Physik, Technische 
Universitat Dresden, D-01062 Dresden, GermanJ*] 

(Dated: March 10, 2015) 

Abstract 

We study the Hilbert-Schmidt measure on the manifold of mixed Gaussian states in multi mode 
continuous variable quantum systems. An analytical expression for the Hilbert-Schmidt volume 
element is derived. Its corresponding probability measure can be used to study typical properties 
of Gaussian states. It turns out that although the manifold of Gaussian states is unbounded, 
an ensemble of Gaussian states distributed according to this measure still has a normalizable 
distribution of symplectic eigenvalues, from which unitarily invariant properties can be obtained. 
By contrast, we find that for an ensemble of one-mode Gaussian states based on the Bures measure 
the corresponding distribution cannot be normalized. As important applications, we determine 
the distribution and the mean value of von Neumann entropy and purity for the Hilbert-Schmidt 
measure. 
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I. INTRODUCTION 


The analysis of spaces of quantum states equipped with some physically motivated mea¬ 
sure has been a lively field of research in the last decade. These investigations help to shed 
light on ” typical “ properties of these spaces and their elements, such as the purity, or the 
entropy of the states. In the case of multi-partite states, an understanding of metric proper¬ 
ties allows for a quantitative characterization of entanglement [1], Detailed knowledge of the 
distribution of a class of states may also serve to specify a priori probabilities for quantum 
state estimation [2] and for quantum tomography Q, and contributes to the understanding 
of the concept of ’’unknown quantum state“ 4]. 


Most results in this field are restricted to quantum systems with a finite dimensional 
Hilbert space, such as many qubit systems, and cannot be easily generalized to systems with 
infinite dimensional Hilbert space. The important continuous variable systems belong to the 
latter and contain the class of mixed Gaussian quantum states, which will be considered in 
this article. Gaussian states are quantum states that have a Gaussian Wigner function, and 


they can be characterized by their covariance matrix a 
for continuous variable quantum information science 


one 


0 . 


These states are fundamental 


For pure states there exists a unique unitarily invariant measure, the Fubini-Study mea¬ 
sure that is induced by the Haar measure on the unitary group. This measure on the pure 
Gaussian states has already been studied in [7]. In our work presented here, the Hilbert- 

Decause it is one of the most 


10]. To this aim a brief in- 


Schmidt measure for mixed Gaussian states will be examined, 
widely used invariant measures for mixed quantum states [8 
troduction to Gaussian quantum states is given in section QTJ The Hilbert-Schmidt volume 
element for Gaussian states is derived in section IIII Al and statistical properties of states 
distributed according to this measure are calculated in section IIII Bl and IIII Cl In particu¬ 
lar, the distribution of symplcctic eigenvalues is calculated. These values can be used to 
determine all unitarily invariant properties of Gaussian states, such as the von Neumann 
entropy and the purity. In the following section IIVI these results are compared to the Bures 
measure for the simple one mode case, which leads to qualitatively different results. Finally, 
in section [V] a summary of the results and an outlook for further research is given. 
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II. GAUSSIAN QUANTUM STATES 


Consider N Fock spaces 'H l in each of which one bosonic creation operator a] is dehned and 
an orthonormal basis is given by the number states {|n) i |n G No}. An TV-mode continuous 
variable system is a quantum system described by the Hilbert space Ti = 1 'Hi, arising 

from the direct product of these Fock spaces. This space is spanned by the product basis 
(K, ■ ■ -,n N ) = \n 1 ) 1 <g> ... <8) \n N ) N \n u ..., n N G N 0 }. 

The characteristic function of a state p in T~L is dehned as: 


N 


y(a) = tr (p-Dq,), D a = J^Jexp (cqa- — a*a*) a G C 


N 


( 1 ) 


i= 1 


re characteristic function is called the Wigner quasiprob- 
111 ] . A Gaussian state in "H is a density operator with a 


The complex Fourier-transfornr of t 
ability distribution in phase space 
Gaussian characteristic function, i.e.: 

x(oi) = exp(i d T ol — -a T Ea), a = (Reai,..., Re« 7 v, Imaq,..., Iiuccat) 7 ^ G R 2N , (2) 

with the displacement vector d G R 2N and the real, symmetric 2 N x 2 N covariance matrix 
E. The displacement d of any Gaussian state can always be brought to zero by a local 
unitary transformation {^] . Therefore it does not contain any information about correlations 
between subsystems or unitarily invariant properties of the state, and will be set to zero in 
the following. A Gaussian state is then exclusively characterized by its covariance matrix 
E. The corresponding density operator of the state is given by [uj: 

/ exp(-ia T Ea)T)_ Q d 27V a. 

Any matrix E is the covariance matrix of a Gaussian state if and only if 


( 3 ) 


12 


14 ]: 


E + i J > 0 , with the synrplectic form J = 


0 1 
-1 0 


( 4 ) 


where > 0 means positive semidehnite. Gaussian unitary transformations Uq are unitary 
operators of the form: 


N 


U G = exp ( - i ^2 PijCLiCLj + P*j<ij<i] + Qija\aj ), P = P T , Q = Q ] . (5) 

*.i =1 
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Together with global phase transformations the Gaussian unitary transformations form the 
class of unitary operators which transform Gaussian states pv; into Gaussian states pg [7]. 
The transformations Uq are a representation of the real symplectic group Sp(2iV, R), and 
lead to a symplectic transformation S of the covariance matrix AH, 

IAgPtU'g = with S T JS = J. 


14]: 


( 6 ) 


For any covariance matrix £ there exists a symplectic transformation S such that S T Y,S is 
diagonal, with each diagonal entry appearing twice fl3| : 

JV = diag(z/ 1 ,...,z/ A r) 


S US' = D = 



ZA- > 1. 


( 7 ) 


The N values zy are called symplectic eigenvalues of £. They characterize a Gaussian state 
up to unitary transformations, and are in this sense the equivalent to the eigenvalues of 
the density matrix in quantum systems with finite-dimensional Hilbert space. The density 
operator corresponding to the diagonal covariance matrix D is a tensor product of thermal 
states of the harmonic oscillator [5j: 


N 


Pd = 


£ 


u, 1 

i =1 * n=0 


Vi 


\n) i (n 


i 5 


( 8 ) 


V{ + 1 . 

n=u 

huJl '. This state is pure if and only if all z/j are equal 


where one may identify i\ = coth , 9ksT 
to one, that is if it is the vacuum state. Hence all pure Gaussian states (the set of squeezed 
states) are Gaussian unitary transformations of the vacuum state. 


III. HILBERT-SCHMIDT MEASURE FOR GAUSSIAN STATES 

A. Hilbert-Schmidt volume element 

In the literature a variety of different physically motivated measures on quantum state 
spaces is used. A thorough introduction to these measures for systems with a finite dimen¬ 
sional Hilbert space can be found in ]lj. One important and widely used measure is the 
Hilbert-Schmidt measure. It is the measure induced by the Hilbert-Schmidt metric. The 
Hilbert-Schmidt distance of two states p and p' is defined as: 

d H s(p,p') = ^/tr((p-p') 2 ) . (9) 
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An important property of this metric is that it is invariant under unitary transformation 
of the density operators. For two Gaussian states ps and px>, the Hilbert-Schmidt distance 
can be calculated (see Appendix [AJ) : 


d-HsipY,-, PE') — 


1 1 
+ 


\ + \/det£' y/det^E + E') 


( 10 ) 


Since the metric is unitarily invariant, this expression should be symplectically invariant on 
the level of covariance matrices. This is indeed the case, because all symplectic matrices S 
have a determinant of one jldij]: 


det {S t TjS) = det S T det £ det S = det £ 


( 11 ) 


The infinitesimal distance element of the Hilbert-Schmidt metric takes the simple form: 

d s 2 hs = tr((dp) 2 ) . (12) 

If one restricts oneself to the manifold of Gaussian states, dp becomes: 
dPs = ^ f d 2Ar aT>_ a d(exp(-ia r £Q:)) 

= ~n J d 2Af aT>- a exp(-iQ: T £Q:)( - ia T d£a) ( 13 ) 

r21V 

If 1 1 

= — / d 2Ar a/l_ Q! exp(--a r £Q:)(exp(--Q: r d£Q:) - l) = p s+ ds - Ps ■ 

7 r JV J Z Z 


The infinitesimal distance element for Gaussian states ds^- 5 = d HS (pY+dY, Ps) 2 is derived in 
Appendix [Bj 


dsl/o — 


(tr(£ _1 d£)) 2 + 2tr((£ _1 d£) 2 ) 


(14) 


HS 16\/det £ 

The volume element of the Hilbert-Schmidt measure can then be obtained by hireling the 
explicit form of the metric tensor g: 


As 2 hs = £ dE„ 9 „^dE^ AV„s = pdetg P[ dE„, 


( 15 ) 


a,(3 


where a and /3 denote some indices labeling all matrix elements of d£. 
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B. Distribution of symplectic eigenvalues 


For quantum states of systems with a finite-dimensional Hilbert space the distribution of 
the eigenvalues of density matrices for the Hilbert-Schmidt measure can be derived Jsj]. In the 
same way it is possible to derive the distribution of symplectic eigenvalues {vj\i — 1,..., N} 
of Gaussian states for this measure. 

Any covariance matrix £ of a Gaussian state can be written as a symplectic transformation 
of a diagonal matrix D (see equation (ITj) ). One may write: 

£ + d£ = S T DS + dZ = S t (D + 6D)S. (16) 

The expression D + 5D can be decomposed as an infinitesimal shift in the symplectic eigen¬ 
values followed by an infinitesimal symplectic transformation: 

D + SD = (1 + d X t )(D + dD)(l + dX) = D + dD + d X T D + DdX . (17) 


X is the generator of a symplectic transformation and is therefore a Hamiltonian matrix 
( JX) T = ,JX 15]. Inserting in (fT4l) gives (see Appendix ICl) : 

1 = f (trfZr'dD)) 2 + 2tr((ZT MD) 2 ) 

(18) 




16\/det D 

+ 4tr((dX) 2 ) +4tr(dXD~ 1 dX T D)) . 


Since no terms of the form dUdX appear, the metric tensor is of block form g = g D ®g s , with 
g D corresponding to the symplectic eigenvalues and g s to the symplectic transformation. 

A general symplectic transformation S G Sp(2X, R) is of dimension N(2N + 1). Together 
with the N symplectic eigenvalues, the expression S T DS = £ is X-fold overdetermined, 
since the dimension of £ is just N(2N + 1). Thus the symplectic transformation S which 
transforms D to £ is not unique. The symplectic transformation: 


S = QS, Q = 


Q i Q 2 

Q 3 Qa , 


Q t Q = 1, Qi diagonal 


(19) 


also transforms D to £ = S T DS, because D and Q commute: 


S t DS = S t Q t DQS = S t Q t QDS = S t DS . (20) 

The transformations S and S are in this sense equivalent. The matrices Q form a group 
that is the direct sum of N rotation groups SO(2), and are therefore of dimension N. 
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Hence they are the only freedom in the choice of S. In order to derive the Hilbert-Schmidt 
volume element explicitly, a 2 N 2 -dimensional parameterization of infinitesimal symplectic 
transformations has to be found, in which this freedom is not included. 


Any Hamiltonian matrix A" can be written as 


W 


X = 


A B 
C -A T 


B = B 1 


C = C 7 


( 21 ) 


And the generator of the matrices Q is: 


R = -R 1 = 


0 G 
-G 0 


G diagonal. 


( 22 ) 


For a given infinitesimal transformation S = 1 + dA", consider the equivalent transformation: 


, 0 dG \ dA dfl 

1 + dX = 1 + d R)S = 1 + dR + dX = 1+ \ + 

—dG 0 / IdG -d A T 


(23) 


Choosing dG = diag(dcn,..., dcjvjv) and assigning d bu = d (bn + SijCu) gives a suitable 
2A r2 -dimensional parameterization of dA": 


dX 


( dA d B 
ydG -dA T 


d B = d B t dG = dG T (dG)** = 0; i = 1,..., N. (24) 


Thus the volume element d Vhs becomes: 

N N 


N N 

d Vhs = \J det g ^ | du^j ^ | da/ m ^ ^ | db 

q>p =1 


(25) 


2=1 


\,m= 1 


r>s =1 


The metric tensor g is derived in Appendix O The result is: 

/ A ± y/N + l ( l Ar(jv+f)—i 
Vdetg= aN2 - J _[_[ (jg - u m ) . 

l>m= 1 


v nL^' 


(26) 


The integral over the symplectic Group is divergent, since the symplectic group is not a 
compact group. Hence the Hilbert-Schmidt volume of the space of mixed Gaussian states 
is infinite. However, the integral over the symplectic group only contributes as a trivial 
factor, since it is independent of the integrand \/det g. Integrating over the symplectic 
eigenvalues gives a finite value for all N, so that the distribution of symplectic eigenvalues 
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Pn(vi, ■ ■ ■, vn) oc \/det g is normalizable. For the simplest cases N — 1 and N = 2 one 
obtains: 


PM = 2 V 


P 2 (i/i,z/ 2 ) = 


2\2 


525 (vf — v%) 
~8 vfv f 


Vi G [1, oo). 


(27) 


These distributions are illustrated in figure CD In general, large values of symplectic eigenval¬ 
ues occur with a small probability (Hilbert-Schmidt measure) - pure states are more likely. 


This is nicely in line with the corresponding findings for qubit-systems (see [l]). Even the 
node structures (e.g. zero along v\ = zq for N = 2), known from qubits, are reflected. For 
N — 1 the maximum of the distribution is located at the pure states v — 1. This changes 
for N > 2, because the distribution is zero when two symplectic eigenvalues are equal. For 
the case N = 2 the maximum is at = 1, ^ 2 (i) = V%- 




(b )N = 2 


FIG. 1. Normalized distribution of symplectic eigenvalues for Hilbert-Schmidt measure for one 
mode (N = 1) and two mode (N = 2) Gaussian states. 


C. Distribution of purity and entropy 

With the help of the distribution (1251) it is possible to calculate expectation values of 
all unitarily invariant properties of iV-mode Gaussian states for the Hilbert-Schmidt mea¬ 
sure, as all these properties can be expressed solely by the symplectic eigenvalues. From 
the derivation of the Hilbert-Schmidt distance of two Gaussian states (Appendix 0 one 
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immediately sees that the purity // =tr(p 2 ) of a Gaussian state is: 


MPs) = tr ((Ps) 2 ) = 


ydetS 


n N 

i=i”i 


(28) 


A second important quantity, the von Neumann entropy S' of a Gaussian state, is given by 

Q: 


S(f>z) = 


N r .11 

-tr(p s hi ps) = ^ | hi 


U + 1 


Z4 — 1 , (Vi — 1 

m 


(29) 


2 / 2 V 2 

For low mode numbers N the mean values of these quantities can easily be calculated 
numerically with high precision using the derived volume element. The results are shown in 
figure El In general, Gaussian states in the Hilbert-Schmidt ensemble become more mixed 
for high mode numbers, since the purity decreases and the von Neumann entropy increases 
as N increases. This is also the case for many-qubit systems, where the mean purity has 
been computed analytically in Hilbert-Schmidt measure as a function of the Hilbert-space 


dimension. The results can be found in 


10](£gure 2(a)), and can be compared to figure [2] 


It is also possible to obtain the entire distribution of the purity of Gaussian states in Hilbert- 
Schmidt measure analytically (see Appendix [E]). For the more complicated von Neumann 
entropy such an expression cannot be given analytically. For N — 1 the purity is simply 
distributed like ^/Ji. This distribution and the ones for N = 2 and N = 3 modes are shown 
in figure [3] The variance of the distributions is decreasing as N increases. It is important 
to note that as a consequence of the infinite Hilbert-space dimension the purity can be zero, 
however the Hilbert-Schmidt ensemble does not contain infinitely mixed states, so that the 
distribution of the purity is zero for // = 0. 

IV. BURES MEASURE FOR SINGLE MODE GAUSSIAN STATES 


The Bures distance of two quantum states can be defined as a function of the fidelity F: 


dii(p,p') = d 2 (l - l/dft /))) F U ’./>') = Hv^Tb) 


(30) 


An explicit formula for the fidelity of two A-mode Gaussian states has been derived in 
For N = 1 this expression simplifies to: 

2 


F{pTnPT,') = 


det(E + Z7) + F — 


V = (dets- l)(det£' - 1). 


(31) 
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FIG. 2. Average purity (a) and average von-Neumann entropy (b) of a Gaussian state (Hilbert- 
Schmidt ensemble) as a function of the number N of modes. 



FIG. 3. Distribution of the purity of Gaussian states in the Hilbert-Schmidt ensemble for mode 
numbers N = 1, 2,3. 


Following the same steps as for the Hilbert-Schmidt measure, one can calculate the infinites¬ 
imal volume element of the Bures measure in the N = 1 case. For N > 1 this is not possible 
offhand, due to the complicated structure of the general formula for the fidelity. 

Using the notation of the last section, the result is: 


dV* = - 


4 v 2 + 1 s/u 2 - 1 


du da d&. 


(32) 


In Bures measure the distribution of symplectic eigenvalues is not normalizable, as it only 
decreases with 1 jv for large v. Thus the Bures ensemble is much more concentrated on mixed 
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states than the Hilbert-Schmidt ensemble. One can define the mean value of a quantity f(y) 
with respect to the not normalizable Bures distribution in the following natural way: 


(f)B 


lim 

I'm ^OO 


JT m mPsM± 

/r PB(v)dv 


PbW) = 


V 




(33) 


It is easy to see that with this definition the mean purity of one-mode Gaussian states in 
the Bures ensemble (fj) B is zero, and the mean von Neumann entropy (S) B infinite. Thus 
when no further restriction on the states is made (e.g. fixed energy), the Bures ensemble of 
Gaussian states contains almost only infinitely mixed states. It is not clear whether this still 
holds for N > 1. However, a consequence might be that the probability to find an entangled 
Gaussian state in the Bures distribution is zero, as numerical calculations suggest for N = 2 

Q. 


V. SUMMARY AND OUTLOOK 

In contrast to systems with finite dimensional Hilbert space, much less is known about 
the geometry of quantum state spaces in continuous variable systems. Gaussian states allow 
for an easy access to this topic because of their simple structure. In this work the Hilbert- 
Schmidt measure has been studied on the domain of mixed Gaussian states. It turns out 
that in spite of the fact that the manifold of Gaussian states is unbounded, Gaussian states 
distributed uniformly according to the Hilbert-Schmidt measure (|25|) have a finite mean von 
Neumann entropy and a nonzero mean purity, which is reflected in the normalizability of 
the distribution of symplectic eigenvalues. This is not obvious, as it turns out that for an 
ensemble of one mode Gaussian states for the Bures measure (1321) this distribution cannot 
be normalized, and the ensemble consists of almost only infinitely mixed states. It is not 
clear whether this is also true for the Bures ensemble of multi mode Gaussian states. 

With the distribution of symplectic eigenvalues at hand, similar to the distribution of eigen¬ 
values in qubit systems, all unitarily invariant properties of the ensemble can be calculated. 
As examples, we study the distribution of purity and the von-Neumann entropy of Gaus¬ 
sian quantum states. Entanglement does not belong to these unitarily invariant properties. 
However, the Hilbert-Schmidt volume element (TT5l) can be used straightforwardly for the 
study of entanglement in multi mode Gaussian quantum states. For any compact subset 
of two mode Gaussian states (for example the setcounterof two-mode states for fixed en- 
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ergy expectation value), the probability for a separable state can be calculated numerically 
for the Hilbert-Schmidt measure using the analytical form of the volume element and the 
Simon-Peres-Horodecki separability criterion jl9|. In addition, it would also be possible 
to consider conditioned bipartite Gaussian states and study ” conditioned volumes”. In this 
area some new results for the two qubit system have been obtained recently [9], and it would 
be interesting to see whether some of those results still hold for Gaussian states. 
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Appendix A: Hilbert-Schmidt metric for Gaussian states 

The expression ([9]) is to be calculated for Gaussian states. 

tr(Os - Ps0 2 ) = tr((p s ) 2 ) + tr((p S /) 2 ) - 2tr(p s piy) 

The characteristic function of the product of two Gaussian states has been calculated in 


17]. The trace is given by evaluating the characteristic function at zero. One obtains: 

1 


tr(psps') = 


det±(E + E') 


from which equation (fTOl) follows immediately. 


Appendix B: Distance element for the Hilbert-Schmidt metric 

The infinitesimal distance element for Gaussian states is: 

—i —i dS _i 

ds HS = dms(ps+d£,Ps) 2 = (dets) 2 + ( det(E + dE)) 2 - 2( det(E + —)) 2 

= (detS)“^l + (det(l + S- 1 dS))“^-2(det(l + E- 1 ^)) - ^ . 

With use of the relation det(M) = exp (tr(ln(M))) 20] and the series representation of the 
matrix logarithm. Up to second order in dE one obtains: 

(det(l + E _1 dE)) 2 = exp ( — -tr(ln(l + E'MS))) 

= 1 - yr(E-ME) + itr^E-ME) 2 ) + i(tr(E-ME)) 2 . 
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The result (TT4|) follows from the same calculation for the second determinant. 

Appendix C: Distance element for the symplectic eigenvalue decomposition 

Using ([TED one obtains for the two terms in (HTTP : 
tr(£ -1 d£) = tr^-MD) + 2tr(dA), 

tr((£ _1 d£) 2 ) = tr((D"MD) 2 ) + 2tr((dX) 2 ) + 2tr(D~ 1 dX T S T DdX) 

+ 4tr(dDdXD - 1 ). 

From J 2 = — 1 and (JdA") T = ,/dX we find: 

tr(dX) = —tr( J 2 dX) = -tr( J( JdX) T ) = -tr( JdX T J T ) = -tr(dX) = 0 . 

In addition [J, D] = [J, D~ x } = [J, d D) = [D, d D) = 0 and D 7 = D applies, and therefore: 

tr(dDdXU)' 1 ) = -ti(dDJ 2 dXD- 1 ) = -tr(dDJ(JdX) r U>- 1 ) 

= -tr(dDJdX T J T D- 1 ) = tr(dD J 2 dXD -1 ) = -tr{dDdXD ~ x ) = 0 . 

Thus all terms proportional to d/MA" vanish. Overall one obtains the result (fTHlh 


Appendix D: Explicit form of the metric tensor 


The metric tensor g D arises from the terms in (TT8T) containing d D. Explicitly, they are 
given by: 


(tr(£)- 1 d£))) 2 = 4f;5^ 2tr((ZT 1 d J D) 2 ) =4^ 

ij =1 * ^ z=l 

Thus the explicit form of g D reads: 

d = 1 (1 + 8ij) _ 1 (1 + 8jj) 

13 4 Vdet D 4nf =1 ^ Wi 

Using the parameterization (124|) . the terms containing dX give: 

N N 

tr((dX) 2 ) = 2tr((dA) 2 ) + 2tr(d5dC') = 2 £ daijdciji + 2 (dbijdcij + dc tJ db t 

i,j =1 i>j =1 


N 


d ± 
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tr(d XD~ 1 dX T D) 

= tr(dAJ\f~ 1 dA T J\f) + tr(dA T Af~ 1 dAN') + tiidBAf^dBAf) + tr(dOA/’ _1 daA0 


N N N 



Overall, the distance element can be written as: 



One can read the explicit form of the metric tensor. The measure y/detg turns out to be: 


m- 1 ( 1 \ N{N+l * ] 

4»d n " =i J 

1 / 1 \N(N+i) 


'«(^) II (£+£)’-< 


* J l>m= 1 ^ 1 


4JV2V nf =1 


If AT 2 ^ ) ( -i—rJV 

n fc=1 u k v m.=i^ 


2(JV—1) 


, 2 _ 7/2 \2 

y Z OrJ 


4JV2V nf=i 


JV(JV+f) —1 


/ ivTT JJ 


.2 , ,2 \ 2 


l>m= 1 
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For any given N the integrals can easily be computed analytically. 
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